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Abstract 

In this paper we are concerned with the effects of anisotropic pressure on the 
boundary conditions of anisotropic Lane-Emden equation and homology theorem. 
Some new exact solutions of this equation are derived. Then some of the theorems 
governing the Newtonian perfect fluid star are extended taking the anisotropic 
pressure into account. 

1 Introduction 

The Newtonian theory of stellar structure is usually used to describe the stars 
with not extremely high densities IB in m Even for high density compact 
objects like white dwarfs and neutron stars, applying Newtonian gravity leads to 
acceptable results comparable to the results of the relativistic models [5]. 

The contents of the star is often modelled by a perfect fluid with an equation 
of state relating the mass density p and pressure P. Among all possible choices of 
this equation, the polytrope equation of state, P = Kp 1 is of great importance for 
many astrophysical situations mm® where K and 7 = 1 + ^ are both constants 
called polytropic constant and polytropic exponent, respectively, (n is called the 
polytropic index). The polytrope equation is the simplest equation of state useful 
for a wide range of fluid densities. Moreover it leads to the scale symmetry and 
therefore introducing homology invariants [4j. 
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On the other hand there are some phenomena leading to the local anisotropy of 
pressure inside a self gravitating system. For low density objects, the Newtonian 
approximation is valid. Some of the physical mechanisms for anisotropy in this 
regime are |6]: 

• When there is an anisotropic velocity distribution in a collisionless gas [7j, 
the radial pressure satisfies Jeans equation: 


dP r 

dr 


d4> 2 

= P~r + -A(r) 

ar r 


a) 


where <fr is the Newtonian gravitational potential of the fluid and we have 
considered a static and spherically symmetric distribution of matter. A = 
P t — P r , where P r and P t are the radial and tangential pressure respectively. 
In this case A measures the anisotropy of velocity distributions: 


A = p((V 2 ) - (V 2 )) (2) 

where (V 2 ) and (V 2 ) are radial and azimuthal velocity dispersions. (Note 
that according to the spherical symmetry assumption we have {V 2 ) = ( Vq )). 
Equation (PQ) describes the hydrostatic equilibrium of a self gravitating ob¬ 
ject. 

To produce anisotropic velocity distribution, consider a galactic halo of 
fermionic dark matter. The conservation of angular momentum of neutrinos 
streaming into the halo, leads to anisotropic pressure [8j. 

• For a slowly rotating system HE], the equation of hydrostatic equilibrium 
in the first order is given by equation m in which A = — |/xu 2 ?’ 2 where oj is 
the angular velocity. 

• For a mixture of two non interacting perfect fluids, the energy-momentum 
tensor is that of an anisotropic fluid in which p, P r and Pt are some functions 
of the mass densities and pressure of each component |9j. 

• For a low mass charged white dwarf, the repulsive electrostatic force can be 
regarded as a source of anisotropy [TU] . 

According to the above cases, even in the Newtonian regime, generally we are 
dealing with two components of pressure. 

Here we shall focus on the Newtonian anisotropic polytropes. In the next 
section, first we derive the anisotropic version of dimensionless Newtonian hydro¬ 
static equation, Lane-Emden equation and it’s boundary conditions. We extend 
the homology theorem for anisotropic polytropes and then in section 3 we de¬ 
rive a new set of analytical solutions. There are some exact solutions of the 
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Lane-Emden equation for special values of barotropic index in arbitrary spatial 
dimension [I]. In this paper, we assume that the form of the Lane-Emden equa¬ 
tion, doesn’t change after introducing the anisotropy factor. Thus we find the 
anisotropic modified form of the existing isotropic solutions by this ansatz. The 
authors of mi. propose another ansatz on A in the modelling of relativistic stars. 
This is used in [I] to generate some anisotropic Newtonian polytrope solutions 
from isotropic solutions. Other various assumptions on A can be found in [2]. In 
section 4, we focus on some integral theorems on the physical quantities of a star 
in the hydrostatic equilibrium obtained by Chandrasekhar in his book pj. These 
are derived from Newtonian equations directly without any special assumption 
on the internal structure of the star or it’s equation of state. These theorems 
contain some inequalities on the central pressure, the gravitational potential and 
the mean gravitational acceleration of a star. Here we extend these theorems to 
the general anisotropic equilibrium configuration. 


2 Anisotropic polytropes 

The hydrostatic equilibrium of a star is governed by equation (JT]) . The Newtonian 
potential is related to the density of the fluid by Poisson’s equation: 

V 2 </> = -4? xGp. (3) 


Considering the fluid obeys the polytropic equation of state along with the equa¬ 
tion (HD, one gets: 


4>(r) -<j)Q = I\(n + 1) 





2A(x) 

xp(x) 


dx 


( 4 ) 


where po and 0o denote the central density and gravitational potential. We take 
<f >o = 0. Substituting equation Q) into equation @ leads to the following expres¬ 
sion which is the fundamental equation of equilibrium: 


K(n + 1 )V 2 p™ — 


1 _ c _[_ ( r N- i2A(r) 

dr y rp(r) 


-N—l 


= —AirGp 


( 5 ) 


where N is the dimension of space. 

Introducing the new dimensionless variable £ 
the following definitions 


± («+l)A^ A 2 

47rGpo 1- " / 


r, and with 


P = P O 0 n , P r = P 0 O n+1 


( 6 ) 
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equation (J5J) takes the following form 


6 " + 


N-K, 


£ 


P 0 {n + 1)&■ 


' , A - 2 A 0' A 
A +—— A-n-A 
s 0 


±0 n . 


( 7 ) 


where prime denotes differentiation with respect to the new radial coordinate, 
£, and the plus and minus signs correspond to —oo < n < — 1 , —l<n< +oo 
respectively. Note that 0 = 0(r) is a dimensionless function. Equation dTj) is the 
anisotropic version of the well-known Lane-Emden equation. 

For n = — 1, the polytrope equation of state gives P r = K if p / 0. Thus ac¬ 
cording to equation ([U) , the term involving A only contributes to the gravitational 
potential and the Poisson’s equation gives: 


A = “T§T J o xN l p( x )dx (8) 

This shows that the case with n = — 1 is a special case which must be treated 
separately. Whereas for the isotropic case (A = 0), one gets p = 0 and therefore, 
n = — 1 is excluded from the further study for isotropic fluids. 

The other important case arises when n = Too and therefore the polytrope 
equation of state reduces to 

P r = Kp. (9) 


Replacing this into (P) , we obtain 


= - r-^dx ( 10 ) 

\PoJ Jo p{x)x 

The combination of the above expression when (f >o = 0 and (J3J) gives 

9" + ^^9' + -^-e e (A' + ^^A + Ad') =e~ e ( 11 ) 

? Pqq Vs / 

where p = p 0 e ~ 9 , P r = P 0 e~ e and £ = ( 4 ^;) 2 r. 

Either the Lane-Emden equation (JTj) (for n / —1, ± 00 ) or (fill) (for n = ± 00 ), 
contains two unknown functions, 0(£) and A(^). Thus we need another equation 
to obtain a closed system of equations. 

Before describing our procedure to introduce the other equation, let us assume 
that A is a given function. To obtain a unique solution for 8 one has to specify 
two boundary conditions. The first is 0(0) = 1 (for n / —l,±oo) and 0(0) = 0 
(for n = ± 00 ). In order to obtain the other condition, let’s integrate (|3|): 


— = —r 1-Ar [ 4:TTGpx N ~ 1 dx ( 12 ) 

dr Jo 
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( 13 ) 


Combination of the above relation with (jT]) and ((6]) , produces 

d0_ __J_i 2 (A(Q\ 

d^° N^ 0+ P 0 (n + l)«ToV £ )■ 

A similar calculation for n = ±oo leads to: 


#L\ = 1 , 2 m 

d£ l^° N^° P 0 \^o £ 


(14) 


Relations m and m lead to the extensions of the Chandrasekhar’s theorem 
[3] for anisotropic polytrope: 

Theorem I: 

The finite solutions of anisotropic Lane-Emden equation at the origin have to 


satisfy f |^=o = 




P 0 (n+1) llm €^0 l 


( *P) 


if n / -1, ±oo and | ?=0 = 4- lim^ 0 




A (0 


if n = ±oo. 

Another important theorem is a generalization of the homology theorem [4]: 

Theorem II: 

If the anisotropic Lane-Emden equation is satisfied by #(£) and A(£) then 

2 2(~+l) 

A™-!0(A£) and A !-« A(A£) satisfy the equation dTJ) , and 0(A£) — In A 2 and 
^■A(A£) satisfy the equation (fTTT) . where A is a constant. 

The proof is straightforward. It can be easily obtained by substituting these 
expressions in (JT]) and (HID . Thus we obtain a set of solutions parametrized by A. 

One usual approach to solve 0 (or (fill) ) with boundary conditions obtained 
in theorem I, is to suppose a special form for A. Then putting it into 0 (or 
mi) yields a differential equation for 6. If this procedure leads to smooth and 
non-negative pressures and density for the fluid, then they can be regarded as the 
physical solutions. However, knowing the function A, equation 0 (or m) is a 
non-linear differential equation which in general one does not expect to obtain an 
analytical solution. Here we introduce a heuristic approach to get some analytical 
solutions. In order to do this, we assume that the anisotropy factor has no 
influence on the form of isotropic Lane-Emden equation. This means that A 
modifies only the coefficients of isotropic Lane-Emden equation. Hereafter, we 
adopt the minus sign on the right-hand side of Lane-Emden equation 0 (or m) 
since the plus sign leads to an imaginary radial coordinate [3]. In the following 
we shall consider two cases: 

Case I: 

A modifies the coefficient of O' in the isotropic Lane-Emden equation . In this 
case equation 0 can be considered as two separated equations as below 


-2 


P 0 (n + l)O r 


' . N — 2 0'' 

A + —— A-n-A 
? 0 


= N\6' 


o" + Nl+N ~ l e' = -t 


(15) 

(16) 
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where Ni is an arbitrary constant. Equation m is a first order differential 
equation with the solution 

(n + l)P 0 iVie n 


A(0 = - 


o, e AT ~ 2 di 


(17) 


2£ n ~ 2 j 

Equations (fl6l) and (1171) with the boundary conditions (fl3l) and 0(0) = 1, give 
the solutions of anisotropic polytropes. 

Case II: 

A modifies the coefficient of 9 n in the right hand side of the isotropic Lane- 
Emden equation. 

-2 


P o (n+l)0 n i 


A/ N - 2 A 0 . 
A + —— A-n-A 

d 2 e n -i(W _ 

dr/ 2 rj dr] 


= N 2 9 r 


where r] = (N 2 + 1) £, N 2 > —1. The solution of equation (fTKl) is 


A(0 = 


~(n + l)P 0 N 2 9 r 


/ ' 




(18) 

(19) 

( 20 ) 


2(JV 2 + l)^" 2 

Now let us to interprete physically the equation (f7|) . The quantities of physical 
interest, for definite values of the polytropic index, are the stellar radius and mass. 
These are: 

i 

/ fiTfr) 4- 1 m 2 1—n 

( 21 ) 


R = 


K(n + 1)Y 

Po Ci 


AnG 


M 


= 1 4 " 2pdr= 41 (^sr 1 ) 1 Pa C reii 


( 22 ) 


where is the first root of 9 function and thus this point defines the surface of 
the star. Using the anisotropic Lane-Emden equation (JTj), we get: 


M = — 47t 


K{n + 1) \ 2 Y d 


4t tG 


Ydf^dd 2 

po it U ~ prr— 

Jo <% V Po{n + l) 


d£ (23) 


which with the boundary conditions 0q = 1, the relation (11311 and the definition 
of surface £i gives: 


M = An 


I\ (n + 1) 2 

AnG 


Po 


,d9 

‘d£ 


(24) 


£i 


This has the known form of the mass relation for isotropic star. Eliminating the 
central mass density between (1211) and (1241) we have: 


3 —n 

MR ™- 1 = 47 t 


K(n + 1) 
AnG 


n+1 

Ci "' 1 


d9_ 


(25) 


£i 
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Figure 1: Plot of as a function of -jj- in the case 1, for n — | with N\ = —1.5 
(dotted line), N\ = 0 (thin line, the isotropic case), N\ = 0.5 (dashed line), N\ = 1.5 
(thick line) and N l = 2.5 (dot-dashed line). 


knowing the polytropic index, 




can be obtained from (1161) or (1191) numerically 


and thus the constant on the right hand side of equation (1251) is found. Therefore 
the anisotropy factor doesn’t appear explicitly in the mass or the mass-radius 
relation of a star. The effect of anisotropy is included in 6 function through it’s 
equation (parameter N\ in (fl6l) and variable ?? in (fT9l) h Figures (P) and (P show 
the nrass-radius diagrams for low density white dwarfs with n = | obtained from 
equations (IIHI) and (fT9l) respectively. 
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Figure 2: Plot of 


in the case II, for n = § with N 2 = 


-M- as a function of ^ , 

Mq Kq 7 2 

(dotted line), N 2 = 0 (thin line, the isotropic case), N 2 = 0.5 (dashed line), N 2 
(thick line) and N 2 = 2.5 (dot-dashed line). 


- 0.8 
= 1.5 










3 Exact analytical solutions of anisotropic Lane- 
Emden equation 


3.1 Case I 

3 . 1.1 ?^ = 0 


With n = 0, the density is constant, p = po, and equation (fTBl) takes the following 
form 

9" + + 1 = 0 (26) 

in which S = N\ + N , where solution is 


9 = 


2S + 1 ’ S > 1 

— 25 + ^ 2—5 "hi) S < 1, S ^ 0 


(27) 


in which A is a constant. In the above solution we have also implemented the 
boundary condition 0 ( 0 ) = 1 and only those terms that lead to the non-singular 
solutions have been written. 

We can now read A from (1171) considering the boundary condition (1131) : 


(S-N)Pp t-2 
2SN ? 


S> 1 


A = 


(S-N)Po 

2 




S < 1,5 / 0,5 + N 


(28) 


0 N=S =1 

We see that, to satisfy the boundary conditions, we have to have A(0) = 0. The 
resulting A function is plotted in Fig.©. Also the radial and tangential pressures 
are determined from the expressions ©, m and ( 1281 ) . 


3 . 1.2 n — 1 


In this case, the differential equation (fl6l) reduces to 

9" + ^ S ~ ^ 9' + 9 = 0 

The solutions of (1291) are given in term of Bessel functions as follows 


m = c 


S — 2 
2 


c\Js^(0 + c 2 Ys^($) 
2 2 


(29) 


(30) 
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Figure 3: Plot of A for S' > 1 in three dimensions with S — 4 (thick line), S — 5 (dotted 
line), S — 6 (dashed line), S — 7 (dot-dashed line). 
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Figure 4: Plot of A for n — 1 in three dimensions with S = 4 (thick line), S = 5 (dotted 
line), S — 6 (dashed line), S — 7 (dot-dashed line). 

where C\ and Ci are integration constants. Imposing the condition 0(0) = I, Ci 
vanishes and we get 

0(O = r (f)(f) 2 ^(0, S > 2 (31) 

in which the constant C\ is expressed in terms of Gamma function, T (^). Sub¬ 
stituting (13TD into m leads to 

m = ^(s-n)p, b« 3 -4Jsja(or 2 (|)r(^)F (f; i + f. i + f i -j ) (32) 

where F(a,b,c,x ) is the hyper geometric function. These results also satisfy the 
boundary condition (fl3l) with 6'{ 0) = A(0) = 0. The function (|32l) is depicted in 
the Fig. dH). 
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3.1.3 n=fr| 

In this case, if one takes 

m = 


n— 1 




(to — l ) 2 

then equation (fTBl) reduces to 

d 2 z z=fz n 
dt 2 (■n — l ) 2 

The above equation can be integrated as follows: 


£ = e 


-t 


1 - 4 
dt) (n — l) 2 


?n +1 


2 T Z 

T T n + 1 


+ C 


(33) 


(34) 


(35) 


Considering n > 1 (S' > 2), turning back to the original variables and using 
equations (l33l) and (l35l) , a simple calculation shows that the initial conditions on 
6 leads to z = 0 and dz/dt = 0 at £ = 0. This gives (7 = 0. Writing the above 
equation in terms of the original variables, we have: 


2 00 ' 6' 2 6 n+1 

7 £“7" + £ TT — ^ 

n — 1 2 to + 1 

We proceed with integrating the above equation 

<9^ 

t— = D 


(36) 


(37) 


The constant D might be determined taking into account (1131) and 0(0) = 1 


1 

15 


_ s_ 6' 

lim 0 s ~ 2 — 

£ 


1 , 2 A(Q 

N + P 0 (n +1)1^0 e 


(38) 


As we see, the constant D depends on the A function near the origin. The 
differential equation (1371) has the following solution 


0 = 



e 

D(S- 2) 


S—2 

2 


(39) 


Employing (1171) . we get the following expression for the difference between the 
radial and tangential pressures 


A(£) = 


SPo(S-N) 2 
ND(S — 2) 5 



£ 2 

D(S- 2) 


(— — l + — -—-^ 

V 2 ’ 2’ 2 ’ D(S-2)J 


(40) 

Now, we are in the position to fix the constant D. Combining (1381) and (1401) . one 
gets D = —S. The resulting A is shown in Fig. ©. 


12 





















Figure 5: Plot of A for n — in three dimensions with S 
(dotted line), S — 6 (dashed line), S — 7 (dot-dashed line). 


= 4 (thick line), S — 5 
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Figure 6: Plot of A for n = 0 in three dimensions with N 2 — 1 (thick line), N 2 = 2 
(dotted line), N 2 = 3 (dashed line), N 2 = 4 (dot-dashed line). 


3.2 Case II 

3 . 2.1 n = 0 

The solutions of differential equations and (l20|) satisfying the boundary con¬ 
ditions are: 

-»2±±e + AVN 2 + lZ + l, N = 1 


e = 


Na+lrf I i 

2N “ r 


A = 


'2{N 2 +1)N 


e, 


N > 1 
N> 2 


2(iV 2 +l) 

The behaviour of A is shown in Fig. 


P n N 2 ( c2 V(^2 + l) 3 ]\T — 1 

2(JV 2 +1)^ N 2 JV - 1 


(41) 


(42) 
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Figure 7: Plot of A for n = 1 in three dimensions with N- 2 — 1 (thick line), N 2 = 2 
(dotted line), N 2 — 3 (dashed line), N 2 = 4 (dot-dashed line). 


3 . 2.2 n = 1 


The situation for this case is also similar to the previous one and the substitution 
of S and £ by N and r/ in (|3Tj) produces the solution of (fl9l) with n = 1. This 
outcome with the integration (1201) give the difference between the radial and 
tangential pressures: 


A(£) = -2 


N-2 


PoN 2 




, ^ (N +^ 2 (-)e- N j^^N 2 +1 ojk(vn 2 + 1 £) 

(iV 2 + 1)2 (7V+1) 2 2 2 


(43) 


which are depicted in Fig. CD- 


3 . 2.3 


n = 


N+2 

N-2 


Similar to the two former cases, in this situation the solution to the (1191) with 

N+2 
N-2 

N 2 + 1 


n = might be obtained by changing S and £ to N and q in (1391) . 


1 - 


D(N — 2) 


N-2 

2 


N > 2 


(44) 
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Figure 8: Plot of A for n — jr^ in three dimensions with N 2 — 1 (thick line), N 2 = 2 
(dotted line), N 2 = 3 (dashed line), N 2 = 4 (dot-dashed line). 

Putting the above result in (1201) gives: 

A _ np 0 n 2 

(N 2 + 1)(N - 2)N 

Combining (|38j) and (l45l) , one gets D = — ■ The graph of (l45l) is is plotted 

in Fig. ® for various values of N 2 . 

4 Integral Theorems 

Chandrasekhar in his famous book [3] has discussed some inequalities for the 
physical quantities describing a star in the Newtonian gravitational equilibrium. 

In this section we will extend some of his results for anisotropic stars. The results 
are general and are not restricted to the polytropic case. 

Let m(r) be the mass contained inside radius r, then 

rn = Anpr 2 (46) 


1 - 


1 + N 2 
D(N - 2) 


e 2 




(45) 
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Equation (pQ) can now be written as 


, Gmp 2 A 

r r ^2 ' ^ 


(47) 


where we have used the definition of gravitational potential. From (|46l) and (1471) 
one can get the following equation: 


1 


r 


2 



—AirGp 


These equations can be used to prove the following theorems: 

Theorem 1: 

For any equilibrium configuration the function 


(48) 


Gm 2 r A 

I = P r + - T — 2 

8 vrr 4 


r, A 

! / dr — 

Jo r 


does not increases outward. 
Proof: Let’s calculate 


Then by (147)) 

QED. 


T . , Gmm! Gm 2 2A 

I = P-\ - 

47 rr 4 27rr 5 r 


T . Gm 2 
I = -f < 0 

27rr 5 


Corollary : 

For the central pressure, we have 


(49) 


(50) 


(51) 


P 0 >Pr + 


Gm 2 
8ir r 4 



> 


GM 2 

8ttFP 



(52) 


in which M = m(R ) is the mass of the star. The last term in the right hand side 
shows a lower bound on Po- Moreover for any arbitrary radius, the above relation 
leads to: 


P r > 


G_ 

87r 




(53) 


Theorem 2: 

For any equilibrium configuration 


if v < 4. 

Proof: 


h 



dm 


Gm 


r R 

4 -7T / dvr 3 " [(4 — v) P r + 2A] 

Jo 


(54) 
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From (1471) we have 


Gmm! 2 A 


4-7T r 4 


- P' 

r 


Multiplying this by r 4 v and then integrating from r to R gives: 


rR 


r -R 


= 8 ir drr 1 ' 6 A — Air / dP, 


~A-u 


(55) 


(56) 


Integrating by parts the first integral, leads to (1541) . QED. 
For v = 4, equation (1561) gives: 

f R A 

I 4 = 8tt / dr— + 47rPo 

Jo 


(57) 


which also has a lower bound equal to ^4 according to 

The gravitational potential energy 14 of the configuration is —I\. According 
to (1561) it is given by: 

f R 

fi = - / dV (3 P r + 2A) (58) 

Jo 

Denoting the mean value of gravitational acceleration by g, we have: 


rR 


Mg = 


dm 


Gm 


= h 


(59) 


and hence by (1561) : 


rR 


rR 


Mg = 8 ir drr (P r + A) = 87 r / dffPt 


(60) 


This shows that for an anisotropic star only the tangential pressure contributes 
in finding the mean value of acceleration. 

Theorem 3: For any equilibrium configuration 

rR 


'kuPqPJ v + ---- + 27 ti /(4 — v) I drr 3 v 


8 R V 


!0 


'0 


R A 

dr — 
r 


f R A GM 2 

> 27 tv I dr n + I u > 


1 0 r 


is —3 


2 By 


( 61 ) 


Proof: 

By theorem 1: 
GM 2 Gm 2 


8 ttR 4: 87rr 4 

and by theorem 2 : 


f R , A f r , A Gm 2 ^ 

— 2 / dr— + 2 / dr— < P r < Pq — ——j + 2 


/o 


/■* 

47r(4 — i/) / drr 3 v 

Jo 


1 0 


GM 2 Gm 2 


87rr 4 


r r A 

/ dr— (62) 
do ^ 


87ri? 4 87rr 4 


[ R A 

f r , A' 

/ dr-h 2 

/ dr — 

Jo r 

/o r . 
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rR 


< I u + 8tt I dff 3 V A < 


This inequality can be written as: 


f R 

47r(4 — is) ( 

Jo 


4ttP 0 R 4 - u + 


rR 


dff 


rR f r 

87tz/(4 — is) dff 3 ~ u / df— 

Jo Jo ? 


3 —v 


Po ~ 


Gmr 
87r r 4 


+ 2 


r df— 

Jo r . 


(63) 


> I v + 8n / dff 6 v A + 


3—y a , I 4 -") _ CM 2 


■"+ 1 > 1R“ 


Inserting 


^ ,.Gm 2 GM 2 
dr — — = —rr-2i„ 


^ —|— 1 


B? 


in ([Ml) and simplifying it, we get (RUT) . QED. 

Corollary : Setting v = 1, we have 


ttPqR + 


3 GM 2 
8 R 


r R 

+ 6ir dff' 


[df± 

Jo r 


(64) 

(65) 


/■ R _2 a ^ GM 2 

>2ir drr 2 A — 11 > 


/o 


2i? 


( 66 ) 


which shows the upper and lower bounds of potential energy, —H. 


5 Concluding remarks 


In this paper we discussed how the anisotropy factor modifies the Lane-Emden 
equation and homology theorem. We obtained some theorems governing the 
characteristic functions of anisotropic star. These are the extension of Chan¬ 
drasekhar’s theorems. We performed a procedure to find the anisotropy factor, 
the radial pressure and the density functions exactly satisfying the Lane-Emden 
equation. We had two cases. In the first case the effect of anisotropy is to change 
the dimmension of space from N to N + N± (which can be a non-integer number). 
For the second case, anisotropy shows itself as a rescaling of the radial coordinate 
in the function 6 . 

Moreover it is straightforward to find some approximate analytical solutions of 
the polytropic stars with anisotropic pressure. For example, let’s assume that the 
including anisotropy factor is equivalent to a slight modification of the polytropic 
index from it’s value in the absence of anisotropy, no- This means that if an exact 
solution do °f isotropic Lane-Emden equation is known: 

0o + ^0o = -C- (67) 


Inclusion of anisotropy factor leads to: 


6 " + 




£ 


P 0 (n 0 + l)£0 n ° 


1 6 r 

A' + -A - uq-^-A 


—9 n ° 


( 68 ) 
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which is equivalent to perturbing 9$ and uq in equation (16711 in the way that 
6 = Oq + £0* and n = no + e satisfy (|67l) : 


N — 1 

9 " + — = ~Q n - 

where |e| -C 1. Thus we have the following equation for A: 


(69) 


2 

P 0 (n 0 +l)^o 


i 6 1 

A' + -A - n o—A 


-9 n + 9 no 


(70) 


Knowing no, the functions 9q and thus 9 can be determined from (1671) and (1691) 
and thus A is obtained from (17011 . With the initial conditions #o(0) = 1 and 
9 q(0), 9*(0), 9[( 0) = 0 and for the case no = 0 and N = 3, we get 9q = 1 — £ 2 /6. 
Inserting 9q into (f69ll , we find that [12]: 


= (3-^)ln(l-^) + ^ln^±i + ^-4 

V 6 ;v 6 J C v/6-e 18 


(71) 


where £ < \/6. Thus: 


A 


Poe 

2 


£ 2 

V[31n(l 


|)-2] + K® arctan _|_4 


(72) 


The plot of A is shown in Fig. ([9]). 
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Figure 9: Plot of A for the analytical solution close to the exact solution (n 0 = 0, 
N = 3) as explained in the text. 
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